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Solutions of the Schrodinger equation are obtained for an electron in a two-dimensional circular 
semiconductor quantum ring in the presence of both external uniform constant magnetic field and 
the Rashba and Dresselhaus spin-orbit interactions of equal strengths. Confinement is simulated by 
a realistic potential well of a finite depth. The dependence of the energy levels on a magnetic field 
strength, strength of spin-orbit interaction and a relative ring width is presented. 
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I. INTRODUCTION 

At present it is well established that motion of an electron in an inner layer of a semiconductor heterostructure 
tIh ' can be treated as a two-dimensional in (x, y) plane because of the existence of the confining quantum well along z 
axis directed perpendicular to (x,y) plane [l|, |2j. In connection with the development of nanotechnology, the study 
of quantum dots and rings in heterostructures acquires increasing importance. In [H, a simple but fairly adequate 
model was proposed in which a two-dimensional circular quantum ring corresponds to an axially symmetric rectangular 
potential well 



!V, 0<p<pi, 
0, pi < p < pa, (1) 
V, Po <p<co 

' of a finite depth V , where p = y 1 x 1 + y 2 , Pi and p Q are inner and outer radii of a ring. Notice, that the confining 
potential of a finite depth was also used for the description of the quantum dots 0, [f| . 

The influence of the spin-orbit Rashba 0, H| and Dresselhaus @ interactions on the electron states in planar 
heterostructures are widely studied in recent years. A uniform constant magnetic field B normal to the plane of 
the quantum ring is described by the vector potential A = y,x, 0). Then the interactions of Rashba Vr and 
Dresselhaus Vd are represented by the formulas 



Vr = a a (a x P y - a y P x )/h, Vd = a D (a x P x - a y P y )/h, (2) 

where P = p + g e A, q e is the absolute value of the electron charge, a x and a y are standard Pauli spin-matrices. The 
strengths of these interactions depend on the used materials. The contribution of two spin-orbit interactions can be 
measured with applying various experimental methods 0, [13] • 

In general, the full spin-orbit interaction has the form Vr + Vd with a R ^ a D . However, the considerable attention is 
spared to the special case 0, [111, [l2[ , when the spin-orbit interactions of Rashba and Dresselhaus have equal strength 
a a = a D — a, which can be experimentally achieved due to the fact that the Rashba interaction strength can be 
controlled by external electric field and the Dresselhaus interaction strength can be varied by changing the width 
of the quantum well along the z axis [H, 0- I n the present paper for the special case a R = a D , we obtain the wave 
functions as well as the dependence of the energy levels on a magnetic field, a spin-orbit interaction strength and a 
relative width of the quantum ring which is described by the potential ([1]) . The calculations are performed for the 
parameter values associated with GaAs. 



* Electronic address: 
t Electronic address: 



kudryash@dragon.bas- net. by 
a.baran@dragon.bas-net.by 



2 



II. EXACT SOLUTION OF THE UNPERTURBED SCHRODINGER EQUATION WITHOUT THE 

ZEEMAN INTERACTION 

The total Hamiltonian of the problem can be written as a sum 

H = H + H' (3) 

of the unperturbed Hamiltonian 

H a = + V c (p) + 2( ffa - a y ){P x + P y ), (4) 

where M e g is the effective electron mass, and the perturbation describing the Zeeman interaction 

H' = \g^Ba z , (5) 

where g is the effective gyromagnetic factor, p, B = My- is the Bohr magneton, M e is the electron mass, <r z is the Pauli 
matrix. 

We shall solve the full Schrodinger equation 

= E^> (6) 
in two stages. First, we shall obtain an exact solution of the unperturbed Schrodinger equation 

H V = E o ^ , (7) 

and then we shall take into account the Zeeman interaction within the framework of the perturbation theory. 

It is easy to see that in the case of the unperturbed equation ([7]) with the Hamiltonian ((4]), in addition to the 
obvious integral of motion 

(8) 



V2 

there is also a non-trivial integral of motion 

L = L z + aM ef t(x-y)(a x -a y )/h, (9) 

where L z is the operator of angular momentum. 

We look for the solutions of Eq. which are eigenfunctions of operators a and L. Then the required solutions 
admit a factorization of the form 

*o(x,2/) = n ± exp \TiV2aM cS (x + y)/h 2 } e im<l> u(p), m = 0, ±1, ±2, . . . , (10) 
where m is the angular momentum quantum number, are eigenvectors of the operator a: 



±( 1 

V2 \±e-™' A 



(11) 



Here we use the polar coordinates p, <f> [x = pcos<^, y = psin0). 
Introducing dimensionless quantities 

p 2M eff p 2 2M cS f% 2M cBPo q eP 2 

r = -, e = ^^E , v = —^V, a = a, b=—B, (12) 
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we write the radial equation 



d 2 u 1 du 
dr 2 r dr 



2 + Z— + ( e a + a 2 ~v c (r) - — -2bm-b 2 r 2 ) u = 0, (13) 



where 



v, < r < ri, 

v c (r) = { 0, n < r < 1, (14) 
v, 1 < r < oo. 

Here we use the notation ri = pi/p for the relative width of the ring. 

In the case of a rectangular potential well (fT4)l in each of the three regions (0<r<r i ,r i <r<l,l<r<oo), 
Eq. (|13[) coincides with the equation for the radial wave function of an electron in a uniform magnetic field without 
taking into account the spin-orbit interaction [l3j . provided we replace e to e + a 2 — v c (r). Therefore, following 13], 
we represent a required function u(r) by the formula 

M f-br 2 \ 

u(r) — (br 2 ) 2 exp f — — j w(r), (15) 

where 

{ciwi(r), < r < n, 

C2iw 2 i(r) + c 22 w 2 2{r), n < r < 1, (16) 
c 3 w 3 (r), 1 < r < oo. 

Here, ci, C21, C22 and C3 are arbitrary coefficients, and the functions u>i(r), it)2i(r), W22(^) and u>3(r) are expressed in 
terms of the confluent hypergeometric functions of the first and second kind M(a, j3, £) and U(a,/3,£) [13] as follows 

Wl (r) =M( 7o ,/3,6r 2 ), 
wai(r) =M( 7l ,/3,6r 2 ), W22W = 17 (7,,/3,&r 2 ) , (17) 
w 3 (r) = U( lo ,/3,br 2 ), 

where 

m + |m| + 1 e + a 2 — i? 
7o = 2 4b ' 

_ m + \m\ + 1 e + a 2 (18) 

(3 = \m\ + l. 

The particular solutions in the first and third regions are chosen so that the radial wave function is regular at the 
origin r — > and tends to zero at the infinity r — >■ 00. 

The continuity conditions for the radial wave function u(r) and its first derivative u'(r) = du(r)/dr at the boundary 
points r — ri and r = 1 lead to a system of algebraic equations 

M/i(m, e , v, a, b, rj)X = (19) 
for four coefficients, where X = {ci, C21, C22, C3} and M^lm, e , v, a, b, r^) is 4 x 4 matrix of the form 



Ma = 



(Wx{n) 


V 



-w 2 i(n) -W22(n) 



-w' 2 i(ri) 
w 2 i(l) 

t"2l(l) 



-w'22 (n) 

w 22 (l) 

^2(1) 






-t«3(l) 
-™ 3 (1)/ 



(20) 
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Therefore, the exact equation for determining e (m,v, a, b, fj) reeds 

det M±{m, e , v, a, b, n) = 0. (21) 

From Eq. (IT3"1) . we see that the dependence of e on a is trivial e (m, o, a, &, = e (m, v, 0, 6, r,) — a 2 . In addition, 
the following relation e (m,u, a,b, r*) — e (— to, v, a, 6, r-j) = 4&to is fulfilled. Of course, Eq. (|2"Tj) cannot be solved 
analytically, but can be easily solved numerically. 

In order to construct the radial wave function completely, we find the values of required coefficients 

-ciMz l {m,e ,v,a,b,n)\ I, (22) 








where 



-w' 21 (ri) -w'^in) 

M 3 = ( Utoi(l) IB22(1) -108(1) | . (23) 
««2l(l) ^22(1) -«/ S (l)> 

The residual arbitrariness in the choice of the coefficient c± is used to implement the standard normalization condition 

(*o ± l*o ± > = l- 

III. CONTRIBUTION OF THE ZEEMAN INTERACTION WITHIN THE FRAMEWORK OF THE 

PERTURBATION THEORY 

The expression for the Zeeman interaction in dimensionless quantities takes the form 

2M oSf % , g M c ff 
* * =4s6<T " * = W (24) 



Since each energy level of the unperturbed system is doubly degenerate with two eigenfunctions (|10j). we consider 
the contribution of Zeeman interaction with the help of the perturbation theory in the degenerate case. 

Because of a z = n T in the basis of the eigenvectors ) and \*S?q) of the unperturbed Hamiltonian, we have the 
following equalities 

(tffHttJ) = 0. (25) 
for the diagonal matrix elements. Off-diagonal matrix elements are given by 

(*+H* ~) = (^-\a z \^)=S(m,v,a,b,n), (26) 

where 



Jo(2ar)u (r)r dr 

5(m, v, a, b, n) = ^ — ^ (27) 

u 2 (r)r dr 



u 



is expressed in terms of the Bessel function. 
Then we get splitting 



e ± = e ± e' (28) 
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for the unperturbed energy levels, where 



e' = e'(m, v, a, b, r^) = 4sbS(m, v, a, b, r.i). (29) 

The following relation e'(— m, v, a, b, r$) = e'(m, V, a, 6, r^) is fulfilled for the corrections e'. Normalized eigenfunctions 
in zero-order approximation, which correspond to the eigenvalues e ± , are described by the formulas 



** = 71 ( * 0+ ± * o) ' (30) 



Note that in the limiting case a = the expressions (28) and ([50)1 become exact. 



IV. NUMERICAL RESULTS 

Now we present the results of the calculations of the quantities e and e'. In accordance with fl5| we choose the 
parameter values M e g = 0.067Af e , g = —0.44 related to GaAs. Then we get s = —0.00737. If we assume p a = 30 nm, 
then the following correspondences a = 1 — > a = 18.9579 meV nm, e = 1 — > E = 0.631933 meV between the 
dimensionless and dimensional quantities are obtained. For example, at chosen parameters, the dimensionless value 
v = 400 corresponds to the potential well depth V — 252.772 meV, which is close to the value 257 meV in Q. 
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FIG. 1: Dependence of e and — e'/e on b at a — 1 and n — 0.5. 




FIG. 2: Dependence of eo and —e/e on a at b = 1 n = 0.5. 

Fig. 1 shows in dimensionless units the dependence of the unperturbed energy e and the relative corrections — e'/e 
to the energy on a magnetic field b for the fixed values a = 1 and r, = 0.5. Fig. 2 demonstrates the dependence of 
e and — e'/e on the spin-orbit interaction strength a for fixed values 6=1 and = 0.5. The solid lines represent 
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TABLE I: Dependence of e and — e'/e on r» at a = 1 and fe = 1. 





m - 


= 


m - 


= 1 


m 


= 2 


e , (-e'/eo) 





4.48334 


26.952 


14.677 


45.9308 


27.3538 


67.4988 




(5.1620 x 10~ 3 ) 


(7.6793 x 10~ 4 ) 


(1.3573 x 10~ 3 ) 


(4.3783 x 10~ 4 ) 


(6.5366 x 10~ 4 ) 


(2.8612 x 10~ 4 ) 


0.1 


7.72294 


36.3199 


14.9000 


47.1566 


27.3612 


67.5757 




(2.7373 x 10~ 3 ) 


(5.5037 x 10~ 4 ) 


(1.3272 x 10~ 3 ) 


(4.2162 x 10~ 4 ) 


(6.5331 x 10~ 4 ) 


(2.8550 x 10~ 4 ) 


0.5 


26.4059 


106.878 


30.3106 


110.949 


37.9733 


119.165 




(5.6869 x 10~ 4 ) 


(1.3905 x 10~ 4 ) 


(4.9374 x 10~ 4 ) 


(1.3401 x 10~ 4 ) 


(3.9012 x 10~ 4 ) 


(1.2497 x 10~ 4 ) 


0.9 


218.124 


404.788 


221.245 


407.21 


226.608 


410.102 




(3.8120 x 10~ 5 ) ( 


-1.2162 x 10~ 5 ) 


(3.7573 x 10~ s ) ( 


-2.0380 x 10~ 5 


(3.6655 x 10~ s ) 


(-1.9963 x 10~ 5 ) 



the first energy levels and the dashed lines represent the second levels for three values of the angular momentum 
(m = 0,l,2). 

Table 1. illustrates the significant dependence of the unperturbed energy e D and the relative correction — e' j e D on 
the relative width r 4 of the quantum ring at the fixed values a = 1 and b = 1. The values of — e'/e are placed in 
parentheses. The value r» = corresponds to a quantum dot. The table lists only first two energy levels. Note that 
for Ti — 0.9 already the second energy levels lie above the potential well with the dimcnsionlcss depth of v — 400 in 
the considered case. In the absence of a magnetic field, the bound states occur only at e„ < v. The presence of a 
magnetic field leads to appearance of the discrete energy levels exceeding magnitude of a potential well depth. 

V. CONCLUSION 

Wave functions and energy levels are obtained for electrons in two-dimensional quantum rings with taking into 
account the spin-orbit Rashba and Dresselhaus interactions of equal strengths in the presence of an external uniform 
constant magnetic field in the framework of an adequate model with a confining potential of finite depth. These 
results may be of interest in the study of spin-dependent effects in semiconductor heterostructures. 
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